Quantum phase transitions in an effective Hamiltonian: fast and slow systems 
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An effective Hamiltonian describing interaction between generic fast and a slow systems is ob- 
tained in the strong interaction limit. The result is applied for studying the effect of quantum phase 
transition as a bifurcation of the ground state of the slow subsystem in the thermodynamic limit. 
Examples as atom-field and atom-atom interactions are analyzed in detail. 
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I. INTRODUCTION 

Frequently, in the process of interaction between two 
quantum systems, only one of them can be detected ex- 
perimentally. In this case, a variety of physical effects 
appear in the process of such interaction which can be de- 
scribed in terms of an effective Hamiltonian correspond- 
ing to the observed system. The simplest example of such 
a situation arises when a fast system interacts with a slow 
system. Then, the fast system can be adiabatically elim- 
inated and the slow system is described by an effective 
Hamiltonian. These considerations were assumed in the 
famous Born-Oppenheimer approximation. A regular ap- 
proach to the quantum dynamics of the observed system 
is provided by the Lie transformation method [l|, |2| ■ The 
advantage of this method consists in the possibility of 
varying the system's parameters, changing relations be- 
tween them, which allows us to describe different physical 
regimes using the same mathematical tool. In particu- 
lar, such an important example as expansion on the res- 
onances in quantum systems not preserving the number 
of excitations can be obtained [![. In this case a generic 
Hamiltonian governing interaction of two subsystems be- 
yond the Rotating Wave Approximation (RWA) can be 
represented as a series in operators describing all possible 
transitions in the system. 

Several interesting features appearing in the process of 
interaction of quantum systems can be realized by study- 
ing evolution of only two generic quantum system with 
one quantum channel. Even in such a simple case we may 
discriminate at least three interesting limits: a) when the 
interaction constant g is much higher than the character- 
istic frequencies of both interacting systems; b) when g is 
smaller than the frequencies of the systems and c) when 
g is higher than the frequency of one system but smaller 
that the frequency of the other one. 

The a) case of very strong coupling should be studied 
carefully, because using the expansion parameter like an 
interaction constant over a characteristic frequency could 
be quite tricky. For instance, the type of the spectrum 
corresponding to the non-perturbed and to the perturbed 



systems can be different: either continuous or discrete. 

The b) case corresponds to a situation where the reso- 
nance expansion is applicable. This particular case leads 
to dispersive-like interactions [4j] . As it was shown in Ref. 
, the evolution is governed by an effective Hamiltonian 
describing a certain resonant interaction and the repre- 
sentation space of the total system can be always divided 
into (almost) invariant subspaces. 

The last case c) possesses a peculiar property: besides 
of finding a corresponding effective Hamiltonian, we can 
also project it out to the lower energy state of the fast 
system, which would never get excited under given re- 
lations between the system's parameters, and thus, de- 
scribe an effective dynamics of the slow system in the 
limit of strong interaction. It is well known that in this 
regime such an interesting effect as Quantum Phase Tran- 
sitions may occur. 

The quantum phase transitions (QPT) are a common 
feature of non-linear quantum systems. Such transitions 
occur at zero temperature and are associated with an 
abrupt change in the ground state structure. QPT are 
related to singularities in the energy spectrum and, at 
the critical points defining QPT, the ground state en- 
ergy is a non-analytic function of the system's parame- 
ters 0. Qualitatively, for a wide class of quantum sys- 
tems, several important properties of QPT can be stud- 
ied in the so-called thermodynamic (semiclassical) limit 
di]. Then, QPT can be analyzed in terms of a classical 
effective potential energy surface 9]. In this language 
QPT are related to the appearance of a new classical 
separatrix when the coupling parameters acquire certain 
values. According to the standard semiclassical quantiza- 
tion scheme and the correspondence principle, the energy 
density is proportional to the classical period of motion, 
diverging on the separatrix, which explains a high density 
of quantum states at the critical points. 

In this article we study effective Hamiltonians describ- 
ing evolution of a generic quantum system X interact- 
ing with a quantum system Y in the case where the 
characteristic frequency of the system X is essentially 
lower than the corresponding frequency of the system 
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Y, ujx i^Y, and the interaction constant g satisfies 
the strong couphng condition: ujx ^ 9 wy. We 
show that, depending on the type of interaction and the 
nature of quantum systems different physical situations 
take place, but generically such effective Hamiltonians 
describe Quantum Phase Transitions in the slow system. 

II. EFFECTIVE HAMILTONIAN 

Let us consider the following generic Hamiltonian de- 
scribing an interaction between two quantum systems: 

H = ujiXo + LO2Y0 + g {X+ + X_) {¥+ + y_) , (1) 

where Xq and Fq are the free Hamiltonians of the X and 
Y systems respectively, and such that wi ^ u)2. The 
above Hamiltonian does not preserve the total excita- 
tion number operator N — Xq + Yq and, in the limit 
liJitL02 3> 5, leads to the appearance of multiphoton- 
type interactions of the form X^YZ" which, under cer- 
tain physical conditions on the frequencies uJi^2, describe 
resonant transitions between energy levels of the whole 
system (see 0\ and references therein). 

The raising-lowering operators X±, Y± describe tran- 
sitions between energy levels of the systems X and Y 
respectively and consequently obey the following com- 
mutation relations: 

[Xo, x±] = ±x±, [Yo, y±] - ±r±. (2) 

We do not impose any condition on the commutators 
between transition operators, which are generally some 
functions of diagonal operators and of some integrals of 
motion [Ni,Xo] = [N2,Yo]^0: 

[X+,X-] - VxoMXo,Ni), (3a) 

- Vy„MYo,N2), (3b) 

where MXo,Ni) = X+X_ and (/)2(yo,^2) = Y+Y^ 
are some polynomials of Xq and Yq respectively (from 
now on we omit the dependence on integrals Ni^2 in 
the arguments) and Vz4>{z) = (j>{z) — (j>{z + 1). The 
objects {Xo,X±) and (Yc, Y+) are known as polynomial 
deformed algebras slpd{2,R) \V\\ . 



Now, we will be interested in the limit where the slow 
system frequency is less than/or of the order of the cou- 
pling constant, wi < g -C 0^2- Following the method 
described in Ref. [3| we can adiabatically remove all 
the terms that contain the fast system's transition op- 
erators, Y±. In particular, the counter-rotating term 
X+Y++X-Y- and the rotating term X+Y- +X-Y+ can 
be eliminated from the Hamiltonian ([T]) by a subsequent 
application of the following Lie- type transformations: 

Ui = exp[e(X+r+ -X_r_)], (4a) 
U2 = exp [e (X+r_ - X+Y^)] , (4b) 

where the small parameters, e and e, are defined by 
e = < 1 and e = < 1. (5) 

LU2 + "-^l ~ 

The transformations (|la|) and (j4bp generate different 
kinds of terms: such as X!j_Y^ + h.c, X'^Y^ + h.c, 
Y^ + /i.e., and X!J -I- h.c. with coefficients depending 
on X{) and Yq- Under the condition a;i, g <^ L02 all the 
rapidly oscillating terms, i.e. those containing powers of 
Y± , can be removed by applying transformations similar 
to ([U, with properly chosen parameters. Then, the ef- 
fective Hamiltonian is diagonal for the operators of the 
Y system. The result can be expressed as a power series 
of the single parameter 5 = g/oJ2 ^ 1- 

It is worth noting that it is not enough that 5 he a 
small parameter for the formal expansion in Q (and the 
subsequent transformations). A balance is necessary be- 
tween the effective dimensions of the subsystems and 6. 
The effective dimensions of the system depend on the 
order of the polynomials (j)i^2, and on the powers of the 
elements X±^q and Kj-n involved in each transformation. 
It was shown before Q, that the powers of the small 
parameters are increasing faster than the powers of X±fi 
and Y±fl, which implies that we can focus on the effective 
dimensions introduced with 

Taking into account the above mentioned considera- 
tions, keeping only terms up to third order in 5 and dis- 
regarding small corrections to the effective transition fre- 
quencies, we arrive at the following effective Hamiltonian: 



Heff = wiXo + c^sFo - 2uji5''V^,-y<^{XQ,Ya + 1) + g5^y<t>y{Yo) {X+ + X-f 

^{gs^^v {MYo)^lMyo - 1)) ix+ + X^t , (6) 

I 

where as 

$(Xo, yo + 1) - Vx„,y„ [01 (Xo) MYo)] , (7) V,„x„^„yJ(Xo, Fq) = /(^o, Y^) - f{Xo + m,Yo + n), 
and the generalized displacement operators are defined for m and n integers. 
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Because the effective Hamiltonian ([6]) is diagonal for 
the operators of the Y fast system, we may project it out 
onto a minimal energy eigenstate of the Y system, \iPo)y, 
substituting Yq by its eigenvalue yo: lo|V'o)y = yo|'0o)y- 

The first order effect then comes from the term ~ 
while the term ^ {X-^-+X^)^ defines a 
fine structure of the effective potential, obtained after 
projecting the effective Hamiltonian ([6]) onto the state 

It is important to stress that, although (5 is a small 
parameter, the effect of the terms ^ (5", n > 1, could 
be in principle comparable with the main diagonal term 
ujiXq, especially if the algebra of X operators describe a 
big subsystem, i.e., large spin or big photon number. In 
this case non-trivial effects such as QPT may occur. 

Now, we may proceed with analysis of the effective 
Hamiltonian ([6]) in the thermodynamic limit, focusing 
on the possible bifurcation of the ground state. 



III. EXAMPLES 

A. Atom-field interaction (Dicke model) 

The Hamiltonian governing the evolution of A sym- 
metrically prepared two-level atoms interacting with a 
single mode of quantized field has the form 

H = ojin + uj2S,_+g{S+ + S-){a'' + a), (8) 

where n = a^a and Sz.± are generators of the {A + 1)- 
dimcnsional representation of the su{2) algebra. 



1. Ejfective field dynamics 

First let us suppose that the atoms form a fast subsys- 
tem so that, 

Xo = n, X+ = a^ X^=a, Yo = Sz, Y± = S±, 

and thus, (l)y{Yo) ~ C2 — Si -\- S z and (j)a:{XQ) = n, where 
C2 = A/2{A/2+l) is the eigenvalue of the Casimir opera- 
tor of the su{2) algebra (integral of motion corresponding 
to the atomic subsystem). 

Projecting the effective Hamiltonian onto the mini- 
mum energy state of the atomic system \0)at, so that 
yo = —A/2, we obtain the following effective Hamilto- 
nian for the field mode: 

Heff = uJin - Ag5 (a -I- a^)^ + gAS^ (a + a^)'^ , (9) 

where wi = cji(l — 2A5^). 

Rewriting ^ in terms of position and momentum op- 
erators, 

Heff = ^{p^ + x^)- 2AgSx^ + AgAS^x^, 



we immediately detect that QPT in this case is re- 
lated to the bifurcation of the effective potential U (x) = 
{uji/2 — 2Agd)x^ + AgAS^x'^ from a single minimum at 
A < 0Ji/{4gS) to a double well structure at ^ > a)i/(4g(5). 
The physical effect associated with this QPT consists of 
a spontaneous generation of photons in the field mode. 
In some sense, the virtual photons, always presented in 
the Dicke model (|S]), are condensed into the real photons 
after crossing the critical point A — Cji/{AgS). It is worth 
noting that this does not happen if the RWA is applied 
to ®. 

2. Ejfective atomic dynamics 

In the opposite case, when the atoms form a slow sub- 
system we have 

Xo^Sz, X±^S±, Yo = h, Y+ = a\ Y- = a. 

Projecting the effective Hamiltonian onto the minimum 
energy state of the field mode |0)/, so that yo = 0, the 
effective Hamiltonian acquires the form 

H,ff = ^iSz - 4.g6Sl + 2iU25^Sl + l&g5^{Sl, Sz}, (10) 

where Coi — uji — 2102^^ — \%gS^ . 

For our analysis it is convenient to perform a 7r/2 rota- 
tion in (jlOp around axis y, transforming the (jlOp Hamil- 
tonian into 

H,ff = -il:^S^~4g8Sl^2uj2&'^Sl-\^g8^{Sl,S^}. (11) 

In the thermodynamical limit we may replace the 
atomic operators by the corresponding classical vectors 
over the two-dimensional sphere, i.e., 

A A A 

Sz — cos 6, Sx — > — sinO cos(f>, Sy — sin(/>sin0, 

and thus rewrite the effective Hamiltonian pT|) as a clas- 
sical Hamiltonian function, 

Hci — — — {uii cos (f> sin 9 2Ag6 cos^ 6 

-Auj2S'^ cos^ (/)sin^ 9 

-^AA^gS^ cos 9 cos (jysin^ 9). (12) 

The first two terms in the above expression describe the 
thermodynamical limit of the Lipkin-Meshkov model Q 
and determine the critical point of QPT, ^ = AAgS/uji = 
1, which again is related to the bifurcation of the ground 
state: a single minimum at sin = 0, cos 9^ = splits 
into two minima at sint/) — 0, cos 6**, — ±-y/l — for 
^ > 1. It is worth noting that the global minimum of 
Hci at ^ < 1 converts into a local maximum for ^ > 1, 
so that Hci{9^,f,) < Hci{9^,). This means that the atoms, 
initially prepared at the minimum of the Hamiltonian 
function, spontaneously change their ground state energy 
at some value of the system's parameters. Classically, 
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this implies appearance of a separatrix, which leads to 
the discontinuity on the energy density spectrum in the 
thermodynamic limit. It is also worth noting that there is 
a loss of the rotational symmetry in this process: the new 
ground state is obviously not invariant under rotations 
around axis x, while the initial ground state is clearly 
invariant under x-rotations. 

It is easy to see that the last two terms in are of 
lower order in the parameter A5 and can be neglected in 
the first approximation for description of QPT at = 1. 

B. Spin-spin interaction 

As a second example let us consider a dipole-dipole like 
interaction, that is, 

Xq = Szi, X± — S±i., Ya = Sz2, Y± = S±2- 

The effective Hamiltonian for the slow spin system (after 
projecting onto the lowest state of the fast spin system 
with eigenvalue —A2/2) takes the form similar to PU)) . 
with 

+16gA2S^St^ + 2AgS^Aj{Sl S,}, (13) 

where a)i = o^i — 2AiuJiS'^ - 20g6^Al. The first two 
terms are dominant for S 1 and describe the Lipkin- 
Meshkov model, so that the critical point is reached at 
^ = AA2A\g5/Cji = 1 in the thermodynamical limit. The 
effect of the rest of the terms in ([T^ is negligible in the 
vicinity of ^ = 1. 

IV. CONCLUSIONS 

We deduce the effective Hamiltonian of a generic slow 
quantum system interacting with another fast oscillat- 



ing system when the total excitation number is not pre- 
served. Analyzing those effective Hamiltonians in the 
thermodynamic limit we have observed a bifurcation of 
the ground state leading to the effect of the Quantum 
Phase Transitions. 

It is interesting to note that, for multidimensional sys- 
tems, when algebraically the X system is a direct sum of 
several non-interacting subsystems, an interesting effect 
of generation of entangled states (in the non-preserving 
excitation case) can be observed. Really, let us suppose 
that in © Xo,± = Xo,±i + Xo,±2, [Xj,i,Xj^2] = 0,j = 
0, ±; then the corresponding effective Hamiltonian (up to 
a first non-trivial order in S) takes the form 

Heff W UJliXos + Xq^2) + i^2Yo + 

gd[{X+^i + X_,i)2 + (X+,2 + X_,2)' 

+ 2(X+,1 + X_,l)(X+,2 + X^^2)]^yCj)y{Y0), 

where we can clearly see that the last term contains 
the operator product ~ X^^iX^_2 which, together with 
quadratic terms in A"-|- 1(2)1 implies a spontaneous gen- 
eration of entangled states of Xi and X2 starting from 
the minimum energy state. This can be corroborated by 
the entangling power measure by considering a uniform 
distribution of the initial factorized states [l2| . Thus we 
can say that, in the regime studied here, entanglement 
can be generated in a bipartite system the vicinity of a 
Phase Transition. 
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